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Plasmonic rod dimers as elementary planar chiral meta-atoms 
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Electromagnetic response of metallic rod dimers is theoretically calculated for arbitrary planar arrangement of rods in 
the dimer. It is shown that dimers without an in-plane symmetry axis exhibit elliptical dichroism and act as "atoms" 
in planar chiral metamaterials. Due to a very simple geometry of the rod dimer, such planar metamaterials are much 
easier in fabrication than conventional split-ring or gammadion-type structures, and lend themselves to a simple analyt- 
ical treatment based on coupled dipole model. Dependencies of metamaterial's directional asymmetry on the dimer' s 
geometry are established analytically and confirmed in numerical simulations. 
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Metamaterials have attracted intense scientific interest in 
recent years for their unusual physical properties rare or ab- 
sent in nature. One example of such properties is giant opti- 
cal activity^ in composite materials containing spiral-like or 
otherwise twisted elements ("meta-atoms"). Recently, planar 
chiral metamaterials (PCMs) were also introduced^. In PCMs, 
the meta-atoms possess two-dimensional (2D) rather than 
three-dimensional (3D) enantiomeric asymmetry (Fig. (T^-b). 
PCMs are distinct from both 3D chiral and Faraday media 
in that their polarization eigenstates are co-rotating elliptical 
rather than counter-rotating elliptical or circular^. This leads 
to exotic polarization properties, e.g., asymmetry in transmis- 
sion for left-handed (LH) vs. right-handed (RH) circularly po- 
larized incident wave without nonreciprocity present in Fara- 
day media. Such exotic properties, combined with the small 
dimensions of planar structures, make PCMs promising for 
polarization sensitive integrated optics applications. 

Several different PCM designs have been proposed so far, 
the most notable example being the asymmetric chiral split 
rings^. Still, this geometry is complicated enough to make 
microscopic theoretical analysis difficult, and it can be shown 
that simpler geometries may be sufficient to achieve the de- 
sired chiral properties. Indeed, simulations reveal (Fig. [T]:) 
that straightening the split-ring segments into two rods (pre- 
serving the segment's length and mutual orientation) results 
in similar transmission asymmetry for LH vs. RH circularly 
polarized light. Rods, unlike more complicated particles, lend 
themselves to a quite straightforward theoretical analysis (see, 
e.g.,^). Hence it is of interest to explore the potential of two- 
rod dimers as meta-atoms for PCMs. In addition, rod-like 
meta-atoms are much easier for fabrication than split-ring or 
gammadion structures, especially in the optical domain. 

In this Letter, we show that double-rod plasmonic dimers 
can function as planar chiral meta-atoms. We find that the 
effective permittivity of a planar rod dimer supports ellipti- 
cal dicroism necessary for PCM effects^. We confirm both 
analytically and numerically that any dimer with distinct 2D 
enantiomers intrinsically exhibits planar chiral behavior, and 
systematically investigate the relations between the strength 
of chiral properties and the dimer' s geometrical parameters. 
By doing so, we propose a simple PCM design that lends it- 
self to easy fabrication and simple analytical treatment. 
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Figure 1 . (Color online) Example of (a) 3D and (b) 2D enantiomeric 
meta-atoms; (c) difference in transmission for LH/RH-polarized 
wave for split-ring vs. "straighened split-ring" (rod dimer) structure. 



We consider the structure shown in Fig. [2^. Since both rods 
in the dimer are elongated, it is reasonable to start with the 
assumption that they respond to the external electric field Eq 
with an induced dipole moment d\ 2 = Si^2Eo where ccj is 
the polarizability tensor of a single rod. For a parallelepiped- 
shaped rod one can follow the depolarization field approach^ 
to obtain polarizability in the form 



(01 



■ajjij^jij (1) 



where the unit vector Jij denotes the rod orientation; fj, coj, 
Yj, and 5^ are geometry-dependent parameters, and stands 
for tensor (dyadic) product. Fig. ^ shows that Eq. ([T} with 
parameters directly obtained from the rod dimensions and ma- 
terial properties provides a decent reproduction for the scatter- 
ing cross-section spectrum of the rod. By letting the param- 
eters fj, (Oj, Yj, and 5^ be adjusted, a nearly perfect best-fit 
between the expression ([T]) and numerical simulations can be 
achieved. Also note that when a rod is bent to represent a 
ring segment, its scattering cross-section near its fundamental 
resonance does not undergo significant changes. 

When there are two rods in a dimer, the external field Eq 
should be modified due to the presence of the neighboring rod. 
This results in a system of coupled equations for the dipole 



dimer considered in the present study, (b) Scattering cross -section 
of a parallelepiped- shaped rod with dimensions 13 x 0.8 x 0.8 mm 
determined by numerical simulations and using Eq. ([TJ, ab initio and 
using best-fit parameters. Numerical results for a 13 mm long ring 
segment are included for comparison. 



moments of rods 

dl,2 Si,2 



Eo + -G(R)d2,i 



(2) 



where R = 7?r is the vector connecting the rods (Fig. |2^) and 
G is the dyadic Green's function in free space: 



G(R): 
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k^R^ J k^R^ 



R^ 



(3) 

with I being the identity tensor. Solving Eqs. (O with (O for 
di 2, one arrives at the total dipole moment d = di +d2 = 
SeffEo of the dimer with effective polarizability tensor 



ai;Ui (8)Jii + a2M2(^M2 + Q^iQ^2?^(Mi (^M2 + M2(K)Mi) 
1 - aia2/c2 

afjii (g) All + C(fji2 (8) Ai2 + 0^3^ (Ml (8) Ai2 + M2 Ml) . 

(4) 

where k defines a coupling coefficient between the rods: 

k^ 

K= — [Gi (Ml • M2) + Gr {Jii • r) {1I2 ■ r)] . (5) 

The effective permittivity tensor Ceff is derived from Sgff in 
Eq. dll) as Ceff = I — (^oVceii)~^Seff. In axial representation, 
Ceff can be expressed as 

^ 1 a^^ 

^eff = I + — TT ^ (C+ C_ + C_ C+) (6) 

^O^^cell ^ 

with the complex vectors c± = (/li + T]±/l2) determining op- 
tical axes and ri± given by 



af ± J {af) - afaf 



a' 



eff 



(7) 



The dielectric permittivity tensor (|6]l corresponds to an ab- 
sorbing nonmagnetic crystal, which in general has two distinct 
eigenmodes with different polarizations, phase velocities and 
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Figure 3. (Color online) (a) Refractive indices for two co-rotating 
elliptically polarized eigenwaves [eigenvectors hi 2 of the tensor Nh 
in Eq. ^] in the effective medium given by Eq. (b) Transmission 
and absorption for incident wave with polarization coincident with 
hi 2. (c) Absolute value of the complex projection coefficients Ci^2 
[see Eq. (|9])] for RH (solid) and LH (dashed) circularly polarized 
incident wave, (d) Analytically calculated difference in transmission 
for LH/RH-polarized incident wave. Here ai = 13 mm, ^2 = 10 mm, 
d = \0 mm, = 45°, and 1// = 0° 



absorption coefficients^. The medium's eigenmodes hi 2 and 
their associated refractive indices ^1^2 are given by eigenvec- 
tors and eigenvalues of the refractive index tensor- 



(- 



-111 



(8) 



Here n is an unit vector of the wave normal, the operator is 
defined as (n^)u = [n x u], and A~ denotes a pseudoinverse 
tensor to A. 

In Fig. [3^ the real and imaginary parts of ^12 are shown 
for propagation direction orthogonal to the meta-atoms plane. 
One can clearly see that absorption bands for the two eigen- 
modes are different, which is characteristic for dichroic me- 
dia. Fig. [3b shows the transmission and absorption spectra of 
the medium given by Eq. Q for an incident wave with polar- 
ization coincident with the eigenmodes (hi ^2). We calculate 
the transmission spectrum from Ceff using generalized transfer 
matrix techniques (see, e.g.,^^). Due to the small variation in 
the real part of the refractive indices, reflection from such an 
effective medium slab is fairly small, and the dips seen in the 
transmission spectra are primarily due to absorption. 

For an arbitrarily polarized incident wave, the slab acts as 
an absorbing polarization filter splitting the incident field H 
into two waves with polarizations parallel to the crystal eigen- 
vectors hi 2, namely 



H = Cihi+C2h2, 



(9) 



where Ci^2 are complex projection coefficients. The transmis- 
sion of the effective medium is then determined by the rela- 
tions (i) between complex projection coefficients Ci^2 and (ii) 
between absorption coefficients for the eigenwaves. 

Fig [3]: shows the projection coefficients Ci^2 for circularly 
polarized incident waves. In this case, it can be seen that 
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Figure 4. (Color online) Numerical (top) and analytical (bottom) dependence of ATrl on (a) inter-rod distance d, (b) difference in rod length 
ai — ^2, (c) rod misalignment angle 0, and (c) rod displacement angle i/A for = 90° (as the dimer changes between T- and L-shaped). Unless 
specified otherwise, ai = 13 mm, a2 = 10 mm, d = 6 mm, (j) = 45°, and y/ = 0°. 



\Ci \ = \C2\. However, in the spectral range of strong absorp- 
tion the coupling of RH and LH incident waves to the crystal 
eigenmodes hi 2 is very different (Cf 2 7^ ^1^2)' leading to pla- 
nar chirality in the form of enantiomeric asymmetry in the 
transmisison ATrl = Tr — Tl (Fig.[3]l). 

Such a difference can be used to quantify the "strength" 
of planar chiral properties in a particular structure. Fig. |4] 
shows the spectra ATrl for various shapes of the dimer, an- 
alytical results from Eqs. O-® compared to the results of 
direct 3D frequency-domain numerical simulationsii. It can 
be seen that the effective medium model offers a good coinci- 
dence with numerical results. 

Most rod dimers are seen to exhibit chiral properties un- 
less the dimer has an in-plane mirror symmetry, which makes 
the 2D enantiomers indistinguishable and therefore enforces 
ATrl = 0. For the geometrical tansformation considered, such 
symmetry is achieved in the following cases: (i) for rods 
of equal length (V-shaped dimer. Fig. lib); (ii) for parallel 
rods (II- shaped dimer. Fig. |4j:); (iii) for a T-shaped dimer 
(Fig. |4]l). In all these cases, the analytical model correctly 
predicts the absence of enantiomeric asymmetry. Otherwise, 
chiral properties appear to be stronger when enantiomers are 
most distinct, so there is an optimum rod misalignment an- 
gle 22.5° in Fig. lit. It is also necessary that the res- 
onances of the individual rods have some degree of spectral 
overlap. Hence, ATrl depends on the length mismatch be- 
tween the rods in a non-monotonic way: it first increases when 
the dimer deviates from the achiral V shape, reaches the max- 
imum value, and then decreases again with a pronounced res- 
onance splitting as Aa becomes greater (see Fig.|4j)). Finally, 
ATrl becomes smaller as the distance between the rods de- 
creases (Fig. 13^). 

Since the coupling between meta-atoms in the 2D lattice 
was neglected both in the analytical model and in the simu- 
lations, it can be concluded that planar chirality in the PCMs 
under study is intrinsic (i.e., attributable to the geometry of 
the meta-atom itself) rather than extrinsic (attributable to the 
meta-atom arrangement, e.g., in tilted-cross arrays ^^). This 
agrees with earlier time-domain simulation results With the 



proposed model extended to include the inter- atom coupling, 
explicit account of intrinsic vs. extrinsic effects in PCMs can 
be given. Combining the intrinsic and extrinsic contribution in 
the same PCM can be used to maximize its chiral properties. 

In conclusion, we have shown that plasmonic rod dimers 
act as elementary planar chiral meta-atoms. Using a simple 
coupled-dipole analytical approach, we show the presence of 
elliptical dichroism and enantiomeric asymmery. We have 
explored the chiral properties of rod dimers of different ge- 
ometries and demonstrated that symmetric dimers are achiral, 
whereas the structures with most distinct 2D enantiomers have 
the most pronounced chiral properties. Analytical results are 
in good agreement with direct frequency-domain numerical 
simulations. 
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